Dark Matter Annihilation around Intermediate Mass Black Holes: an update 
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The formation and evolution of Black Holes inevitably affects the distribution of dark and baryonic 
matter in the neighborhood of the Black Hole. These effects may be particularly relevant around 
Supermassive and Intermediate Mass Black Holes (IMBHs) , the formation of which can lead to large 
Dark Matter overdensities, called spikes and mini-spikes respectively. Despite being larger and more 
dense, spikes evolve at the very centers of galactic halos, in regions where numerous dynamical effects 
tend to destroy them. Mini-spikes may be more likely to survive, and they have been proposed as 
worthwhile targets for indirect Dark Matter searches. We review here the formation scenarios and 
the prospects for detection of mini-spikes, and we present new estimates for the abundances of mini- 
spikes to illustrate the sensitivity of such predictions to cosmological parameters and uncertainties 
regarding the astrophysics of Black Hole formation at high redshift. We also connect the IMBHs 
scenario to the recent measurements of cosmic-ray electron and positron spectra by the PAMELA, 
ATIC, H.E.S.S., and Fermi collaborations. 

PACS numbers: 



I. INTRODUCTION 

The distribution of matter around a Black Hole (BH) 
is inevitably affected by its formation and evolution, both 
in the case of baryons [J, Q and Dark Matter (DM) @]. 
The physics is relatively simple: when a compact object 
forms, the surrounding matter reacts to the increased 
gravitational potential, and if the growth proceeds on a 
timcscalc much longer than the dynamical time of the 
system, large overdensities can be achieved. This is par- 
ticularly important for the DM distribution, because the 
DM annihilation rate is proportional to the square of 
the number density of DM particles. This means that if 
BH formation increases the DM density locally, it also 
boosts the annihilation rate and improves the prospects 
for detecting DM through the observation of secondary 
particles such as gamma-rays, anti-matter and neutrinos. 

Gondolo and Silk have applied this argument to the 
distribution of DM at the center of the Galaxy Q , where 
a Supermassive BH (SMBH) of 2-4 x 10 6 M o is known 
to dominate the gravitational potential within a sphere 
of radius rj, ~ 1 pc. The consequent enhancement in 
the DM density around the central SMBH (dubbed a 
spike) was subsequently shown likely to be reduced by 
major merger events, off-center formation of the seed BH, 
gravitational scattering off stars and DM annihilations [3, 
HI, H ■ Zhao and Silk have subsequently suggested that 
mild overdensities of DM could also be present around 
Intermediate Mass Black Holes (IMBHs), remnants of 
Pop III stars, and provided the first rough estimates of 
the properties of mini-spikes surrounding IMBHs in the 
Milky Way 0. 

Bertone, Zentner and Silk (hereafter BZS) built a con- 
sistent formation and evolution scenario for IMBHs and 
associated mini-spikes, for two different classes of mod- 
els: "mild" mini-spikes around IMBHs remnants of Pop 



III stars, and "strong" mini-spikes, around more massive 
compact objects @. For both scenarios, BZS tracked 
the merger history of each individual BH, and selected 
precisely those IMBHs which never experienced merg- 
ers, to ensure at a minimum that major mergers have 
not destroyed any spike that may have existed around 
the original BH. This has allowed a detailed quantitative 
description of the mini-spikes scenario, and practically 
all subsequent studies published in the literature have 
made use directly or indirectly of the BZS IMBH mod- 
els. Moreover, BZS were able to predict some of the 
demographics of wandering IMBHs, that may not be di- 
rectly associated with the Milky Way bulge or disk, and 
they proposed to search for a class of gamma-rays objects 
sharing the same energy spectrum as a smoking-gun sig- 
nature of DM annihilations. 

Here, we review the formation scenarios and the 
physics behind the BZS IMBH models, and explore how 
the population of IMBHs, thus of mini-spikes, depends on 
cosmological, astrophysical and particle physics parame- 
ters. To this aim, we update the original BZS catalogs, 
including the most recent determination of cosmological 
parameters, and explore the dependence on astrophysical 
quantities. 

We also review here the experimental strategies that 
may allow the detection of mini-spikes, including the ob- 
servation of multiple gamma-ray and neutrino sources in 
the Milky Way a collection of gamma-ray sources around 
M31, and the possible contribution to the absolute flux 
and the angular power spectrum of the Extra-galactic 
Gamma-ray Background (EGB). 

The paper is organized as follows: in Section[Hl we re- 
view the IMBHs formation scenarios and discuss the cur- 
rent population of IMBHs in the Milky Way. In Section 
IIIII we discuss adiabatic contraction and the formation of 
mini-spikes. Section [TVl is dedicated to the prospects for 
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detecting DM through gamma-rays and neutrinos from 
IMBHs in our Galaxy and in M31. In Section [V] we 
discuss the impact of the mini-spikes population on the 
Extra-galactic Gamma-ray Background. In the final sec- 
tion ( Section IVI[) we discuss our results and present our 
conclusions. 



II. INTERMEDIATE MASS BLACK HOLES 
A. Hints for detection of IMBHs 

We refer to BHs with masses within a relatively large 
range, 100M Q < M < 10 6 M Q as Intermediate Mass 
Black Holes. The lower limit is roughly the mass of the 
largest BH that can be produced by stellar collapse [Til ] . 
and the upper limit corresponds to the minimum mass 
for a SMBH. There is not a clear experimental evidence 
for this class of objects, but they are often invoked as an 
explanation of Ultra Luminous X-ray sources, i.e. sources 
that emit at a luminosity larger than 10 39 erg s _1 , thus 
exceeding the Eddington luminosity of stellar BHs 1 1211. 

IMBHs may also be hosted inside globular clusters [l3| . 
A hint in favor of this possibility comes from the fact 
that the predicted masses for the IMBHs and the velocity 
dispersion typical of a globular cluster follow exactly the 
extrapolation at lower values of the M, — a relation valid 
for SMBHs [H[. 

From a theoretical point of view, IMBHs may pro- 
vide massive enough seeds for the growth of high-redshift 
SMBHs, as observed by the Sloan Digital Sky Survey 
in the form of quasars at redshifts as high as z ~ 6 
flU [l6l . Conclusive evidence for IMBHs may come in 
the near future from the detection (e.g. by LISA of 
the gravitational waves emitted in the mergers of IMBHs 

B. Formation 

In BZS, two formation scenarios for IMBHs have been 
studied, and we refer to that paper for a detailed discus- 
sion. In the first scenario (scenario A) IMBHs are rem- 
nants of the collapse of Population III stars [12]. Since 
the Jeans mass scales with the temperature as T 3 / 2 , the 
mass of Pop III stars (which form early in a hot environ- 
ment) is expected to be higher than other stars, around 
150 — 200 Mq [H, [24|. Newtonian simulations suggest 
that the evolution of Pop III stars is very different from 
stars in the local Universe. Zero metallicity Pop III stars 
evolve on a timescale of order t* ~ 1 — 10 Myr. If 
their masses are in the range M ~ 60 — 140Af Q and 
M > 260Mq, they collapse directly to BHs, while if 
140 < M/M G < 260M Q , they are completely disrupted 
due to the pulsation-pair- production instabilit y, l eaving 
behind no remnant [ll] (see also Refs. [H, 13, 11 HI). 

The second formation scenario (scenario B) considered 
by BZS was based on Ref . [3fj| . as an example of models 



in which the collapse of primordial gas in early-forming 
halos leads to the "direct" formation of very massive ob- 
jects [3ll, [H, [H HH, IH, H6|. In this specific scenario, 
the baryons in the low tail of the angular momentum 
distribution form a pressure-supported disc at the center 
of the halo. Gravitational instabilities then turn on an 
effective viscosity that make baryons lose their angular 
momentum, causing an inward mass flow. The process 
comes to an end after 1 — 30 Myr, when the first stars 
form or when the halo experiences a major merger [30l ] . 
The central object will then undergo gravitational col- 
lapse producing the final BH. 

The requirement that loss of angular momentum via 
viscosity is effective, fixes a mass scale of 10 7 — 10 8 Af Q 
for the initial halos, and leads to a mass scale for the final 
IMBHs of order 10 5 M Q . In fact, the mass distribution 
for scenario B turns out to be a log-normal Gaussian 
(a, = 0.9) with a mean value M # : 




where M V i r is the virial mass of the halo, Zf the rcdshift 
of formation, / the fraction of the baryonic mass of the 
halo that goes into the central pressure-supported disc, 
k the fraction of the mass of the disc that ends up in the 
final IMBH and t is the timescale available for angular 
momentum transfer due to the effective viscosity. 

BZS have estimated the population of IMBHs in the 
MW in scenario A as follows: they populated all the halos 
at z = 18 that constitute a 3cr peak (with respect to the 
average density) with a pristine IMBH of mass 100 Mq. 
The evolution of the halos is then simulated as in Refs. 
[l9l . [37l . [H, HI] , a procedure that tracks the growth and 
mergers of the structures down to rcdshift zero. The 
IMBHs in the halos that never experience mergers sur- 
vive until the present epoch and constitute the current 
population of scenario A IMBHs. Their mass is assumed 
to remain unchanged with respect to the initial value 
~ 100 Mq, and the number of unmcrged IMBHs was 
found to be N A = (1027 ± 84). 

In the case of scenario B, BZS have followed a similar 
approach for the evolution of IMBHs, with the difference 
that only those halos with a mass larger than the mass 
threshold M t hr for the onset of the effective viscosity (see 
the discussion above) are populated with IMBHs, with 
a mass function as in Eq. [1] (see also Fig. [T}. The 
predicted number of IMBHs in the MW for scenario B 
was N B = (101±22) i (see below for more recent results 
obtained with updated cosmological parameters). 

Each of these models of high-redshift seed BH forma- 
tion is reasonable, but both of these scenarios remain 
uncertain in their detail. In particular, these models ex- 
hibit parameters that may reasonably take on a wide 
range of values leading to markedly different predictions 
for IMBH demographics at z = 0. The formation of seed 



3 




1 10 100 

Distance to Galactic Center [Kpc] 



FIG. 1: Left: Mass function of unmerged IMBHs in the scenario B for a MW halo at z = 0. The distribution is based on the 
average of 200 Monte Carlo realizations of a halo of virial mass M mr = 10 12 A hT x MQ, roughly the size of the halo of the MW. 
Taken from Ref. [jj. Right: Radial distribution of the IMBH population in the MW from the numerical simulation of Ref. 

The points refers to the average among 200 Monte Carlo realizations of the MW halo and the error bars show the scatter 
among realizations. The solid line is the analytical fit and the dotted line is a NFW profile. Taken from Ref. (TjJ. 



BHs in both scenarios is cut off at rcdshifts when sig- 
nificant small-scale fragmentation of baryonic disks sets 
in and BH formation in scenario B terminates absolutely 
when the Universe becomes significantly reionizcd. In 
what follows we exhibit the significant parameter freedom 
within these models by varying the formation redshift in 
scenario A, and the reionization redshift in scenario B 
among viable and illustrative values. In scenario A, low- 
ering the formation redshift requires that the seed black 
holes generally form in halos of higher mass, leading to 
reduced abundance, while in scenario B lowering the red- 
shift of reionization allows angular momentum loss to be 
effective for a longer period, increasing the number of 
seed black holes. We note here that the range of param- 
eters we choose are currently viable and refer the reader 
to BZS for modeling details. 



The radial distribution of IMBHs was found to be 
steeper than a Navarro-Frenk- White (NFW) profile 
(see the right panel of Fig. [J), and it is well fitted by the 
following analytical function [1 01 ] : 



N(r) = 5.96 x 10~ 2 



9.1 kpc 



-10.8 



kpc 



-3 



(2) 



III. ADIABATIC CONTRACTION 

As mentioned in the introduction, IMBHs and SMBHs 
are particularly interesting for indirect DM searches, be- 
cause they may act as annihilation boosters, as their 
growth can lead to large DM overdensities that enhance 
annihilation fluxes (see e.g. Ref. [4l[ for a review). In 
particular, if the crossing timescale t cr (the time needed 
for a DM particle to cross the central part of the halo) 
is much smaller than the BH growth timescale t growt h, 
DM particles conserve their angular momenta and their 
radial actions (the integral of the radial velocity over one 
closed orbit), and move to orbits which are closer to the 
BH, thus enhancing their number density and annihila- 
tion rate. 

The first study of adiabatic contraction was performed 
by Young in 1980 [|J in the context of stellar systems. In 
subsequent studies also DM has been taken into account 
[42l . l43l and two different situations have usually been 
considered: 

• "analytical cores" , in which case the initial DM 
profile can be Taylor-expanded around the BH as: 

• 7-models, in which the inner initial DM profile is a 
power-law, p x (r) oc r~ 7 . 

The result of the adiabatic contraction is different in 
the two cases. A strong density enhancement, or spike, is 
produced in both cases, but the final slope is steeper for 
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FIG. 2: Possible models for the DM distribution in the Galaxy. The thin curve shows the standard halo model, and the thick 
curve is the same model after "adiabatic compression" by the Galactic baryons (stars and gas). Both curves are normalized to 
a DM density of 0.3 GeV cm -3 at the Solar circle. Curves labeled "spike" show the increase in density that would result from 
growth of a Galactic SMBH at a fixed location. The annihilation plateau, p = p a = m/ (cv)t, was computed assuming m = 200 
GeV, (av) = 10~ 28 cm 3 s _1 and t = 10 10 yr. Dotted vertical lines indicate the SMBH Schwarzschild radius (r Sc h ~ 2.9 x 10~ 7 
pc, assuming a mass of 3.0 x 1O 6 M0 and zero rotation), the SMBH gravitational influence radius (rh ~ 1.7 pc), and the radius 
of the Solar circle (Rq ~ 8.0 kpc). Effects of the dynamical processes (scattering of DM off stars, loss of DM into the SMBH, 
etc.) are excluded from this plot; these processes wouldgenerally act to decrease the DM density below what is shown here, 
particularly in the models with a "spike.". Taken from [7j. 



7-models. The difference is due to the different behav- 
ior of the phase-space distribution function f(E,t). In 
the case of 7-models, / diverges in the limit E $(0), 
meaning that cold orbits (those orbits populated by par- 
ticles with very low velocities) have a large occupation 
number. The final DM profile for a 7-modcl is a new 
power-law p x oc r _7sp with a steeper slope. The value 
of jsp depends on the initial slope as in the following 
equation: 



The effect of SMBH growth at the Galactic center on 
DM density profiles is shown in Fig. [3]for an initial NFW 
profile (7 = 1). The final profile in this case is a power- 
law of index j sp = 7/3 within a region of radius r sp ~ 
0.2771, where rh is the radius that encloses a DM mass 
equal to twice the mass of the central BH. 

There is an upper limit to the DM density that can 
be achieved around the SMBH (see Fig. [2]). In fact, the 
power-law solution is valid only down to a radius equal to 
the last stable orbit around the SMBH and, furthermore, 
DM annihilation itself sets an upper limit on the DM 
density of order m x /crv(t — tf), so that the density is 
saturated at a cut-off radius of 

r cut = Max [4i? Schw , riim] , (4) 

where i?schw is the Schwarzschild radius of the IMBH 
i?Schw = 2.95km(Af bh /M ). 



Any spike forming at the center of a galactic halo would 
inevitably be affected by dynamical processes that tend 
to deplete the DM density jlHll- 

For instance, a merger 
between two halos, the interaction with a globular cluster 
or the presence of a bar are all able to disrupt the dis- 
tribution of the cold particles that constitute the spike, 
efficiently damping the enhancement [4, 5] . Furthermore, 
the BH may not form exactly at the center of the DM dis- 
tribution [J] , in which case the orbits that the BH crosses 
in its spiraling in will be depleted and the final enhance- 
ment will be much shallower than depicted in Fig. [5] 
Finally, one should also take into account the interac- 
tions of the DM particles with baryons Q: transferring 
energy to the DM particles, the effect is, as before, of 
heating the DM so that the particles in the spike will 
be able to leave the central region. When also includ- 
ing the effect of DM annihilations, it is possible to show 
that any DM overdensity around a BH at the center of a 
galactic halo can hardly survive until today @ , although 
interactions of baryons can actually lead to the re-growth 
of mild overdensities known as DM crests [45[ . Our ap- 
proach is to account for some of these mitigating factors 
by considering only those BHs that formed at high red- 
shift and were found never to interact with any other BH 
during the hierarchical merging process, as discussed in 
Ref. Q. Moreover, at least for scenario B, off-center for- 
mation is not possible for IMBHs and interactions of DM 
with baryons have only minor effects, so that, in contrast 
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to SMBHs, IMBHs can effectively retain their mini-spike. 

IV. DETECTION OF DM AROUND IMBHS 

The flux of secondary particles produced by DM anni- 
hilation around an IMBH can be written as 



— (E) 

dE K 1 Hnml dE 



av) dN 1 



(E) 



2 
Psp 



(r)dV , 



(5) 



where (av) is the thermally-averaged annihilation cross 
section, m x is the mass of the DM particle and dN 7 /dE 
is the number of secondary particles produced per anni- 
hilation. The integral extends from the cut radius r cut 
(i.e. the radius where the DM density reaches the annihi- 
lation plateau in Fig. [2|) to the spike radius r sp . Inserting 
the appropriate reference values, Eq. [5] can be recast as 



d$> 

— (E) 

dE y ' 



$0 

dE 
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kpc 

Tent 



(av) 



io- 



,3 B -1 



Psp(j"sp) 



10 2 GcV cm" 3 

5/3 



100 GeV 

2 , v 14/3 
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(6) 



10" 3 pc 



with $ = 9xl0~ 10 c 
makes it evident that IMBHs can be bright DM annihi- 
lation sources, possibly as bright as the entire MW. 

The fact that the cut radius itself depends on the 
mass and the cross section of the DM candidate al- 
ters the naive dependence of the annihilation flux from 
these quantities, and it is easy to show that in this case 
d^/dE oc (av) 2 ^ 7 m^ 9 ^ 7 . The weak dependence of the 
predicted annihilation flux on the particle physics param- 
eters is one of the intriguing features of the mini-spikes 
scenario. 

In Figs. [3] and 0] the number of IMBHs associated with 
a gamma-ray flux (integrated above 1 GeV) larger than 
$ is plotted as a function of <!> itself. With respect to 
the similar Figs. 4 and 5 in Ref. @, we consider here 
the most recent cosmological parameters from WMAP5, 
and find a significant decrease of the number of IMBHs . 
In order to show the effect of changing the cosmological 
parameters we have considered two representative values 
of the redshift of formation Zf for the scenario A and we 
have varied the reionization redshift z r for the scenario 
B in the WMAP5 3a range [4f|. The total number of 
IMBHs for the different choices of astrophysical parame- 
ters is shown in Table fl] 

Figs. [3] and 0] can also be loosely interpreted as the 
number of objects that can be detected with an ex- 
periment of sensitivity For EGRET, $ egret ~ 
3 x 10~ 8 cm~ 2 s _1 , while for Fermi $Fermi ~ 3 x 
10 _10 cm _2 s — 1 . The fact that the scenario is within 
the reach of EGRET, means that some of the uniden- 
tified EGRET sources could be interpreted as IMBHs. 



Scenario 


z 


Nbh 


la scatter 


A 


18 


319 


49 


A 


15 


54 


13 


B 


7 


278 


37 


B 


12 


62 


14 


B 


17 


8.4 


3.7 



TABLE I: Number of unmerged IMBHs in a Milky-Way like 
halo at z = 0. The four columns indicate the formation sce- 
nario, the formation (reionization) redshift in the case of sce- 
nario A (B), the number of IMBHs and the lcrscatter among 
realizations 



If this interpretation was confirmed by Fermi, which is 
expected to see many more sources, the identification 
of a population of point-like sources with the same en- 
ergy spectrum would provide strong evidence for a DM 
annihilation signal. We stress that the number of de- 
tectable sources depend on the formation scenario but 
also on other poorly constrained processes, such as the 
merger history of IMBHs in the galaxy and the initial 
distribution and dynamical evolution of the surround- 
ing DM. Furthermore, from a particle physics point of 
view, we recall that the DM annihilation cross section 
can largely differ from the thermal value in presence of 
efficient coannihilations |47ll . significant Sommcrfeld en- 
hancements [H, |4t| [5^, l5l| . in cases where DM particles 
are not thermal relics, or in non-standard cosmological 
scenarios [H, HI, [H| • 

Even in the conservative case where none of the 
EGRET sources is associated with a mini-spike, there are 
models where the Fermi satellite may detect this class of 
objects, do to the steep decline of the luminosity function 
in Figs. [3jand[2 

The first experimental constraints on the mini-spikes 
scenario have been obtained through a survey of the in- 
ner Galactic plane at photon energies above 100 GeV 
performed by the H.E.S.S. array of Cherenkov telescopes 
from 2004 to 2007 About 400 hours of data have 

been accumulated in the region between -30 and +60 
degrees in Galactic longitude, and between -3 and +3 
degrees in Galactic latitude, and a H.E.S.S. sensitivity 
map was computed for DM annihilations. The data ex- 
clude scenario B at a 90% confidence level for DM par- 
ticles with velocity-weighted annihilation cross section 
av > 10 28 cm 3 s -1 and mass between 800 GeV and 10 
TeV. 

As for the prospects for detecting mini-spikes with the 
Fermi telescope, a dedicated analysis of the minimum de- 
tectable annihilation flux has been presented in Ref. (56j 
and further discussed in the pre-launch estimates of the 
Fermi sensitivity to DM annihilation signals (57j- The 
mock IMBH catalogs obtained for the Milky Way have 
been adapted, by a suitable rescaling, to the population 
of IMBHs hosted by the Andromeda galaxy [58|, lead- 
ing to the prediction that the M31 should host about 
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FIG. 3: IMBHs integrated luminosity function, i.e. number of IMBHs producing a gamma-ray flux larger than a given flux $ 
(above I GeV), as a function of $, for scenario A. For each curve we also show the la scatter among the different realization 
of the Milky Way-size host DM halo. The blue (red) curve refers to a DM particle with a mass of 100 GeV (1 TeV) and a cross 
section of (av) — 3 x 10 _26 cm _3 s~ 1 ((av) = 3 x 10 -29 cm -3 s -1 ) . The upper (lower) set of lines corresponds to a redshift of 
formation Zf = 18 (15). The figure can be interpreted as the number of IMBHs detectable by experiments with a point source 
sensitivity <E> (above 1 GeV) as a function of $. 



65 IMBHs (again this estimate strongly depends on the 
assumed cosmological parameters). Detectability toward 
Andromeda is reduced greatly due to the comparably 
large distance to the Andromeda galaxy. On the other 
hand, in the case of detecting a handful of such objects, 
they will all be located within ~ 3° from the center of 
Andromeda, which may be an interesting spatial signa- 
ture (see Fig. [5]). 



charged current interactions of neutrinos with the nuclei 
around the detector. Accounting for oscillations we can 
write the muon flux as: 



'"■(E) - S / Ar)dV 



dE 



(7) 
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A. Neutrinos from IMBHs 

DM annihilation can also produce neutrinos, both di- 
rectly or through the decay of secondary particles. Be- 
ing strong annihilation " boosters" , IMBHs may also be 
bright neutrino sources [59|,[6(|, which can be detected by 
high energy neutrino telescopes such as ANTARES and 
IceCube or the future Km3. Neutrino telescopes detect 
neutrinos through the observation of muons produced by 



where V(vg — > v^) is the probability of having a 
muon neutrino at the detector when at the source a neu- 
trino of the I family is produced, and it can be writ- 
ten in terms of the oscillation matrix U: V(v£ — > v^) = 

I!.- <..,,.-■. \ u (-i?\ U M?- 

Finally, the rate of events in a neutrino telescope is 



x P v " HN 
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FIG. 4: IMBHs integrated luminosity function, i.e. number of IMBHs producing a gamma-ray flux larger than a given flux $ 
(above 1 GeV), as a function of $, for scenario B. For each curve we also show the la scatter among the different realization 
of the Milky Way-size host DM halo. The blue (red) curve refers to a DM particle with a mass of 100 GeV (1 TeV) and a 
cross section of (av) = 3 x 10 _26 cm _3 s _1 ({av} = 3x 10 -29 cm -3 s -1 ) . The middle (upper, lower) set of lines corresponds to a 
redshift of reionization z r = 12 (7, 17). The figure can be interpreted as the number of IMBHs detectable by experiments with 
a point source sensitivity $ (above 1 GeV) as a function of $. 



where y v = 1 — E™ ,r /E u , E 1 ^ is the experimen- 
tal threshold for muons, A(E fl ) is the effective area 
and P^(E Ul y) is the probability that a muon neutrino 
interacts with a nucleon producing a muon of energy 
Ep = (1 — y)E v (see Ref. I59f for details). Results are 
summarized in Fig. 2 of Ref. [5{j], similar to Fig. 2]in the 
case of gamma-rays. Those results can be easily extended 
to the new models presented here, the main consequence 
being a smaller number of detectable sources. 



V. CONTRIBUTION OF IMBHS TO THE EGB 
AND THEIR ANGULAR SPECTRUM 

As we have seen in the previous sections, IMBHs are 
generically predicted in scenarios that seek to explain 
the existence of SMBHs at the centers of galaxies. It is 
therefore possible, in principle, that the gamma-ray flux 
from DM annihilations around IMBHs in all DM halos, 



and at all redshift, add together and contribute signif- 
icantly to the diffuse Extra-galactic Gamma-ray Back- 
ground (EGB) discovered by EGRET 

The origin of this gamma-ray emission is actually un- 
known, and although unresolved blazars are often consid- 
ered as the most likely sources of the EGB, the most re- 
cent analysis suggest that they can not account for more 
than 25-50% of the measured EGB [6^] . Besides standard 
astrophysical sources, like galaxies or clusters of galaxies 
[H, HE HII , DM annihilations have been extensively stud- 
ied as possible candidates to explain the EGB emission 

The contribution from cosmological DM halos to the 
EGB is constrained to be rather low, in order to sa tisfy 
the observational bounds at the center of our Galaxy [7(| , 
even when the effect of spikes around SMBHs is taken 
into account. The presence of cosmological DM clumps 
could enhance the signal but according to the most recent 
results from DM simulations the expected boost factors 
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FIG. 5: Map of the gamma-ray flux in units of photons 
cm _2 s _1 , from DM annihilations around IMBHs in the An- 
dromeda galaxy M31, relative to one random realization of 
IMBHs in M31. The size of the bins is 0.1° and the thresh- 
old is 4 GeV as appropriate the best angular resolution of 
Fermi (still labeled as GLAST in the figure). The circles 
highlight IMBHs within the reach of Fermi for a 5cr detection 
in 2 months. The big circle shows for comparison the M31 
scale radius r s . 



are small (see Refs. [7l|, uM and also the results in Ref. 
[691]). A sizeable diffuse emission is instead produced by 
the population of unresolved DM clumps hosted in our 
Galaxy (6^, [z3| . 

As for IMBHs, they can account for a large fraction of 
the measured EGB [61( ■ The mean intensity as a function 
of the energy {1(E) dm) produced by DM annihilations 
in cosmological IMBHs is given by: 



(I(E)dm) 



cdz 



W(E[l + z],z), 



(9) 



where 



W(E 



z) = p)dN^ {E[1 + z]) e -T{E { l +z] , z )tf {z y 



8wml dE 



(10) 

and the optical depth r parametrizes the gamma-ray 
attenuation due to interactions with the extra-galactic 
background light (see e.g. Ref. [74|). 

The factor A 2 (z) in Eq. fTOl contains the information 
on the DM spike profile p sp around each IMBH and the 
comoving number density n(z): 



A 2 (z) = n(z) 



P%(r)d 3 r. 



(11) 



The latter input cannot be computed in a straight- 
forward manner using our current simulation techniques, 
but it can be modeled in an approximate manner based 
upon their results. For example, the BH formation mod- 
els assume that initial BH formation is peaked at a cer- 
tain redshift, Zf, to a good approximation, and is effective 
only in DM halos with masses above a certain threshold 
M t hr- These two parameters depend upon the specific 
formation scenario. Following this prescription and em- 
ploying the halo mass function dn/dM(M, z) predicted 
from extended Press-Schechter theory [75[ , the comoving 
number density at the formation redshift is: 



n(zf) 



Mth, 



(12) 



where it is conservatively assumed that only one IMBH 
is formed per halo. 

The number of IMBHs in present halos depend on the 
formation history of the host halo and is affected by the 
occurrence of BHs mergers. In good approximation it 
follows a linear dependence on the mass of the host halo, 
with a normalization N^h for the Milky Way that is ob- 
tained from simulations. The present IMBHs comoving 
number density is thus: 



n(0) 



dM^(M,z 



0)N bh 



M 



M 



MW 



(13) 



For a generic redshift z the IMBHs comoving num- 
ber density can be parameterized assuming a power- 
law behavior which interpolates the BH number den- 
sity at the formation redshift and at z = 0: n(z) = 
n(z f )[(l + z)/(l + z f )f. 

With this recipe and based on the results of the simu- 
lation of Ref. [3| the authors in Ref. [6l[ computed the 
contribution of IMBHs to the EGB for the two IMBHs 
formation scenarios presented in Sec. [ill arid for differ- 
ent choices of DM particle physics properties. In Fig. [6j 
we show that in the conservative case of scenario A the 
expected signal is largely below the EGB measurements. 
However, for scenario B the predicted diffuse emission is 
at the level of EGRET observations already for a stan- 
dard particle physics setup, e.g. assuming a "thermal" 
annihilations cross section (av) = 3 x 10 _26 cm 3 s _1 and 
m x = 100 GeV. 

As with the number of Galactic IMBHs, these predic- 
tions are subject to relatively large astrophysical uncer- 
tainties. Nonetheless, it is remarkable that models with 
astrophysical and particle physics parameters that have 
been fixed for independent reasons, the DM signal from 
cosmological IMBHs is comparable with the EGB flux 
extracted from EGRET observations. Of course, in or- 
der to detect effectively DM annihilations in the EGB it 
is necessary to distinguish the DM signal from other as- 
trophysical emissions, such as that originated by blazars. 
Unfortunately, this is a difficult program to pursue using 
only spectral information, unless the DM signal is accom- 
panied by striking spectral features, like those induced by 
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FIG. 6: Left: Contribution to the EGB from DM annihilations around IMBHs for scenario A and B. m x = 100 GeV, (<tv) = 
3 x 10 _26 cm 3 s _1 and DM annihilations into bb have been considered. Taken from [6lj |. Right: Angular power spectrum of 
the EGB from DM annihilations around IMBHs at an energy Eo = 10 GeV. The dashed line shows the DM contribution, the 
dotted one is for unresolved blazars and the dot-dashed is the cross-correlation term. The black solid line is the total signal 
and the error bars are for 2-years of Fermi data. The solid blue is the signal for DM annihilations in subhalos. Taken from jlQI ] - 



gamma-ray line emission (see Ref. [6l| for computations 
in IMBHs scenarios). 

This problem can however be overcome by studying 
the angular correlations of the diffuse extra-galactic emis- 
sion, as in Ref. [76j |. In fact, the angular power spectrum 
of anisotropies from DM annihilations and from blazars 
arc different, due to their different spatial distribution, 
energy spectra and radial emissivity profiles. Later on, 
other papers computed the angular spectrum of the EGB 
in the case that a substantial amount of the average EGB 
flux is accounted by DM annihilations in extra-galactic 
halos and subhalos [77| or in Galactic DM clumps [73|, 
or considering both the contributions 69]. In all cases, if 
DM plays a role in the EGB, the analysis of the angular 
power spectrum of gamma-ray anisotropies can nicely de- 
tect its contribution and with the soon-available full-sky 
data by Fermi this approach will be testable. 

In the case of IMBHs, the authors of Ref. [lj| focused 
on scenario B and computed the angular power spectrum 
induced by DM annihilations around IMBHs. Wc refer 
to the original reference for details on the computations. 
As shown in Fig. [5] the shape of the angular power spec- 
trum should help in distinguishing the DM signal from 
the blazar contribution to the EGB. These results refer 
to an energy of observations of Eo = 10 GeV, which is 
a good compromise between maximization of the pho- 
ton count and minimization of the Galactic foreground, 
which tends to masquerade the extra-galactic emission. 
In addition, natural DM candidates of around 100 GeV 
produce the largest gamma-ray yield at these energies. In 



order to assess more precisely the prospects for DM de- 
tection in the EGB they considered the unresolved blazar 
contribution as a known background, motivated by the 
fact that it should be modeled quite accurately from the 
future Fermi catalog of detected blazars. In addition, 
the power spectrum from unresolved blazars should be 
energy independent, due to the power-law energy depen- 
dence of their energy spectra and therefore it could be 
calibrated at low energies, where the contribution of DM 
annihilations is negligible and then subtracted from the 
total anisotropy data. After modeling the blazars contri- 
bution to the EGB, if it is assumed that the remaining 
fraction fuM of the EGB flux at the energy of observa- 
tion Eq is due to DM annihilations, the signal Cf depends 
on the DM power spectrum Cf M and on the cross cor- 
relation between the DM and blazars signals Cf r : 

Q = !luG? U + 2f DM (l - f DM )Cf r . (14) 

Fig. [B]shows the signal as well as the projected la error 
bars after two years of Fermi observations. The prospects 
for distinguishing the signal of DM annihilation around 
IMBHs in the angular EGB power spectrum are encour- 
aging. This conclusion holds for different choices of DM 
mass and annihilation channel and energy of detection, 
provided that DM annihilations substantially contribute 
to the EGB at Eq, i.e. for /dm > 0.3. Moreover, if the 
DM signal is dominated by DM emission in cosmological 
subhalos instead that around IMBHs, the predicted an- 
gular power spectrum is different and for sizeable values 
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Bringmann, Lavalle & Salati (2009) 
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FIG. 7: The solid line show the constraint in the DM annihi- 
lation rate F = l/2(av)/m^ J dxp sp (x) of a nearby DM over- 
density at a distance d from the Earth. The shaded area shows 
the combinations of V — d needed to explain the PAMEFA 
positron excess. The dashed lines refers to a set of propa- 
gation parameters defined in Ref. [79(. For comparison the 
dotted line indicates the 10% of F for the whole MW assuming 
(av) = 3 x KT^cm-V 1 . Taken from [§3|- 



of Jdm the two scenarios should be distinguishable by 
Fermi observations. 



VI. DISCUSSION AND CONCLUSIONS 

Indirect DM searches have received a lot of attention 
recently, in large part due to the recent results of exper- 
imental collaborations such as Fermi, PAMELA, HESS 
and others. However, unless unambiguous spectral fea- 
tures are detected [78[ , it will be difficult to discriminate 
any incident flux of high-energy particles from astrophys- 
ical sources (either known or unknown) and identify a 
signal due to DM annihilation. The mini-spikes scenario 
discussed here may provide an interesting solution to this 
problem. 

In fact, in the previous sections we have shown that 
IMBHs are promising targets for DM indirect detection, 
as a consequence of the large DM overdensities produced 
around them by gravitational processes during the for- 
mation of the BH. The properties of these objects and 
the distribution of DM about these objects both de- 
pend on the unknown BH formation mechanism. In Ref. 
[9( the authors investigated in detail two astrophysical 
models proposed in the contemporary literature. This 



study utilized numerical simulations to make predictions 
for the population of IMBHs that should be present in 
our Galaxy today. These objects could appear as bright 
gamma-ray emitters as a consequence of the DM anni- 
hilations occurring in the dense DM spikes. Indeed the 
prospects for detecting such objects with Fermi and Air 
Cherenkov Telescopes are encouraging. A spectacular 
signature for a DM signal may come from the detection 
of a class of gamma-ray point sources that are not associ- 
ated with astrophysical sources and that share the same 
energy spectrum. IMBHs may also be detected through 
neutrinos or they can produce a sizable excess in anti- 
matter cosmic-ray fluxes, as studied in Ref. [sT|] . More- 
over, analogous populations of IMBHs should be present 
in other galaxies, such as the nearby Andromeda galaxy 
and several gamma-ray sources could be detectable. 

In general, the gamma-ray fluxes from IMBHs in all 
DM halos at all redshifts produce a diffuse gamma-ray 
background that for the most optimistic scenario B can 
explain a large fraction of the observed EGB. This opens 
the interesting possibility to detect DM annihilations in 
the extra-galactic background. The best strategy would 
be to analyze not only the EGB energy spectrum but 
also the angular correlations of the emission. In fact, 
the EGB angular power spectrum could provide a robust 
method to distinguish the DM signal from conventional 
astrophysical emission. 

In Ref. [8(|, the authors derived constraints on the 
mini-spike scenario from a comparison of predictions 
based on BZS models with EGRET data. These authors 
claim that they rule out the possibility to detect the prod- 
ucts of DM annihilation in mini-spikes. Despite the bold 
claim, their analysis only rules out a specific combination 
of cosmological and astrophysical parameters, for only 
one (out of many) possible formation scenarios. Actually, 
we showed above that by simply adopting the WMAP 5- 
year cosmological parameters, the expected number of 
objects went down by a factor of 2, and that changing 
the redshift of formation in scenario A, and of ionization 
in scenario B, substantially modifies the predicted num- 
ber of sources. Moreover, Figs. [3] and [4] show that one 
can obtain quite optimistic prospects for Fermi (with up 
to ~ 200 sources in reach of detection) even assuming no 
detection by EGRET. 

Still, the constraints of Ref. [8(| are extremely useful. 
In particular, these results bear on the possible implica- 
tions of the mini-spikes scenario for the recent measure- 
ments of the positron and electron flux at high energy 
[82l . |83| . The possibility that mini-spikes boost the anni- 
hilation signal thus providing the appropriate normaliza- 
tion of the positron flux is found to be viable, as shown in 
Fig. [3 and it would be interesting to see whether such a 
result holds for updated mini-spikes models. This would 
allow to circumvent the multi-messenger constraints (see, 
e.g., Refs. gHH). 

In conclusion, the mini-spike scenario provides the 
opportunity to discover a new population of sources 
(IMBHs) and DM annihilations at the same time. The 



11 



most prominent predicted signature of this scenario is un- 
mistakable: the observation of many point-like sources, 
not correlated with the Galactic disk, with an identical 
spectrum. It is therefore important, despite the various 



constraints and the many astrophysical uncertainties, to 
search for this population of objects in present and up- 
coming data, especially those of the gamma-ray satellite 
Fermi. 
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